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ABSTRACT 


Shock  wave  propagation  and  focusing  theories  were  in¬ 
vestigated  in  the  U  S  Army  Weapons  Command  Science  and 
Technology  Laboratory  to  determine  the  feasibility  of  util¬ 
izing  shock  waves  for  military  applications  The  propaga- 
tion  of  hypothetical  weak  shock  waves  with  small  planar  or 
concave  frontal  areas  was  analyzed  on  the  basis  of  Whitham's 
shock  wave  theory  The  results  indicate  that  shock  waves 
with  these  characteristics  cannot  be  propagated  intact  over 
militarily  significant  distances  The  general  applicability 
of  Whitham  s  theory  to  this  problem  is  questionable,,  however, 
in  view  of  experimental  data  in  the  literature  which  indi¬ 
cate  that  very  weak  shock  waves  can  be  effectively  focused. 
Therefore,  these  discrepancies  should  be  experimentally  re¬ 
solved  to  explicitly  define  the  feasibility  of  military 
applications  of  shock  waves. 
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OBJECTIVE 


The  objective  was  to  investigate  shock  wave  propagation 
and  focusing  theories  to  determine  the  feasibility  of  util¬ 
izing  shock  waves  for  military  applications. 

INTRODUCTION 


For  certain  military  applications,  it  would  be  desirable 
to  focus  shock  waves  or  shocked  sound-waves  into  a  finite 
area  at  a  significant  distance  from  a  source  such  as  a  de¬ 
tonation  tube.  Consequently,  this  investigation  was  initiated 
to  study  the  feasibility  of  accomplishing  a  concentration  of 
shock-wave  strength  or  of  eliminating  the  normal  geometric 
dilution  of  shock-wave  strength  associated  with  shock-wave 
propagation. 

An  approximate  theory  of  shock-wave  dynamics  developed 
by  Whitham*  has  been  applied  to  predictions  of  strong  shock- 
wave  diffraction  and  stability  problems  with  general  success. 
For  shocks  with  low  Mach  number,  however,  the  theory  does 
not  appear  to  be  quite  as  satisfactory  in  the  prediction  of 
observed  shock-wave  diffraction  behavior.2  Nevertheless,  the 
Whitham  theory  appears  to  be  generally  accepted  as  an  ana¬ 
lytical  working  model  for  the  prediction  of  dynamic  shock- 
wave  behavior. 

Recent  data  obtained  in  studies  of  N-wave  focusing  have 
indicated  that  weak  shocks  may  behave  according  to  the  prin¬ 
ciples  of  geometrical  acoustics.3  This  point  would  seem  to 
be  quite  significant  since  the  focusing  behavior  predicted 
by  geometrical  acoustics  is  markedly  different  from  that  pre¬ 
dicted  by  the  Whitham  theory.  These  differences  will  be  con¬ 
sidered  later  in  the  report.  However,  the  shock-wave  pressure 
amplitudes  used  in  the  N-wave  study  were  pertinent  to  sonic 
boom  phenomena,  and  the  largest  amplitude  used  was  of  the 
order  of  one  one-hundredth  of  the  amplitude  of  interest  here. 
This  point  also  would  seem  to  be  quite  significant  since  the 
realms  of  classical  acoustics  and  nonlinear  shock  phenomena 
are  separated  largely  on  the  basis  of  pressure  amplitude. 

The  feasibility  of  shock-wave  focusing  was  primarily 
considered  in  terms  of  Whitham' s  theory  since  it  is  believed 
that  this  theory  should  best  describe  the  behavior  of  shocks 
having  peak  overpressures  in  excess  of  one  psi.  Some  impli¬ 
cations  of  the  geometrical  acoustic  theory  on  shock  focusing 
are  also  made  later  in  the  report,  however.  Focusing  in  the 
usual  acoustic  sense,  in  which  rays  converge  at  a  point,  is 
not  necessarily  required  for  many  military  applications. 
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All  that  is  required  is  to  concentrate  the  major  portion  of 
energy  produced  by  a  generator  into  a  given  area  at  some 
distance  away.  This  could  conceivably  be  accomplished  by 
properly  shaping  the  wave  front  as  it  emerges  from  the  gen¬ 
erator  or  by  some  focusing  device  external  to  the  generator 
The  acoustical  power  output  of  the  generator  is  generally 
distributed  over  some  area  (as  opposed  to  a  point  source) . 
Consequently,  if  it  is  desired  to  concentrate  this  energy  in 
an  area  of  a  few  square  meters  at  a  significant  distance  from 
the  generator  with  minimal  spreading  losses,  the  shape  of  the 
wave  front  should  be  essentially  planar  or  possibly  concave. 
The  general  approach  taken  in  this  investigation  was  that  of 
analyzing  the  propagation  of  planar  and  concave  shock-wave 
fronts  to  provide  insight  into  establishing  the  feasibility 
of  directing  shock  waves. 

THEORY ,  DISCUSSION.  AND  RESULTS 

I.  Planar  Shock  Waves 

A.  Formation  of  Planar  Shock  Waves 

Planar  shock  waves  are  formed  in  a  detonation  tube 
or  in  a  shock  tube  because  shock  waves  propagating  down  a 
channel  become  oriented  normal  to  the  channel  walls.  A  dif¬ 
fraction  process  takes  place  upon  exit  of  the  shock  wave  from 
the  end  of  the  tube,  however,  and  the  details  of  this  dif¬ 
fraction  depend  on  the  geometry  of  the  tube  exit. 

Diffraction  of  shock  waves  has  been  investigated  by 
Whitham  and  by  Skews.  1*2  Predictions  of  the  shape  of  planar 
shock  waves  diffracting  around  convex  corners  were  made  by 
Skews  based  on  the  approximate  theory  of  Whitham.  The  re¬ 
sults  of  these  predictions  were  compared  with  experimental 
results  for  corner  angles  varying  from  15°  to  165°.  The 
theoretical  and  the  experimental  results  were  generally  in 
good  agreement,  and  the  typical  two-dimensional  behavior 
observed  is  shown  in  Figure  1.  In  general,  a  disturbance  on 
the  initial  planar  shock  originating  at  the  corner  spreads 
out  on  the  planar  shock.  The  final  shock  configuration  con¬ 
sists  of  a  short  planar  region  which  is  normal  to  the  new 
wall  direction,  an  expanding  curved  section,  and  the  remain¬ 
ing  portions  of  the  original  planar  shock.  The  angle  m  , 
in  Figure  1,  represents  the  angle  taken  by  the  point  on°the 
incident  planar  shock  where  the  curvature  of  the  disturbance 
begins.  Skews  found  that  experimentally  the  angle  m  was 
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FIGURE  1  SHOCK  DIFFRACTION  PATTERN  AROUND  A  CONVEX  CORNER 


given  satisfactorily  by  the  derived  relation 


t*"2"0  ■  CY-l)  (M*-l)  (M»  •  2/  (vl))  /(ytl)  * 


where  M  is  the  incident  shock  Mach  number  and  y  is  the 
specific  heat  ratio.  A  graph  of  this  relation  is  shown  in 
Figure  2,  with  mQ  shown  in  degrees  as  a  function  of  MQ. 

It  is  of  interest  to  note  in  Figure  2  that  except 
for  very  small  values  of  M  ,  m  is  in  excess  of  20°.  For 
this  reason,  a  planar  shocx-wave  formed  in  a  tubelike  device 
cannot  be  transferred  to  the  atmosphere  intact.  Instead,  as 
the  shock  emerges  from  the  tube  exit,  diffraction  occurs,  with 
the  curved  disturbance  shock  propagating  along  the  incident 
shock  at  an  angle  m  from  the  tube  exit,  as  illustrated  in 
Figure  1„  Since  m  “ranges  from  20  to  30  degrees  for  incident 
shock  Mach  numbers°from  1.1  to  5.0,  the  incident  planar  shock 
can  be  expected  to  be  severely  degraded  at  points  more  than 
about  R  cot  20°  distant  from  the  tube  exit,  where  R  is  the 
tube  radius. 

The  question  now  arises  whether  the  diffraction 
losses  could  be  lessened  if  the  detonation  chamber  exit  geom¬ 
etry  were  somehow  modified.  As  was  mentioned  earlier.  Skews 
studied  the  shape  of  diffracting  shock  waves  for  exterior 
tube  angles  in  the  range  of  15*  to  165°.  The  diffracted  wave 
shapes  remained  essentially  constant  in  time,  and  all  were 
of  the  general  form  shown  in  Figure  1.  Thus,  it  is  not 
apparent  that  any  modification  of  the  tube  exit  angle  can 
significantly  reduce  the  degradation  of  the  planar  wave  shape 
due  to  diffraction. 

The  effects  on  the  emerging  shock  diffraction  process 
of  the  interior  angle  of  the  detonation  chamber  might  also  be 
questioned.  However,  if  a  planar  shock  enters  the  final  sec¬ 
tion  of  the  detonation  tube  traveling  along  the  tube  axis,  a 
converging  or  diverging  final  section  would  be  expected  to 
only  strengthen  or  weaken  the  shock  respectively,  prior  to  its 
emergence  from  the  section  into  the  atmosphere.  When  the 
shock  emerges,  one  would  again  expect  the  diffraction  effect 
appropriate  to  the  strength  of  the  shock  upon  exit  from  the 
section.  The  preceding  discussion  indicates  the  reason  why 
a  pure  planar  shock  wave  of  finite  area  cannot  be  formed  in 
the  atmosphere  from  such  a  shock  traveling  down  a  tube.  If 
such  a  wave  can  be  formed  by  any  means  whatever,  it  would 
apparently  have  to  be  done  by  some  process  external  to  the 
tube  which  is  not  subject  to  diffraction  effects,  such  as 
propagating  a  shock  wave  through  a  nonuniform  medium. 


4 


FIGURE  2  pL0T  OF  n»  =  Tan'1  [(y-lHM^-l)  (M* 


B.  Shock-Wave  Refraction 


A  fundamental  distinction  of  shock  waves  from  sound 
waves  is  that  in  a  given  medium,  sound  waves  propagate  at 
the  prevailing  sonic  velocity,  whereas  the  propagation 
velocity  of  shock  waves  is  dependent  on  the  shock  strength. 
However,  for  a  given  shock  strength  and  a  given  medium,  shock 
velocity  is  directly  related  to  the  sonic  velocity.  Conse¬ 
quently,  in  a  region  in  which  the  properties  of  the  medium 
change  in  such  a  way  as  to  change  the  sonic  velocity,  the 
shock  wave  velocity  changes  correspondingly,  and  refraction 
of  either  sound  or  shock  waves  can  take  place. 

It  is  well  established  that  under  proper  conditions, 
the  atmosphere  can  cause  refraction  of  acoustic  or  shock 
waves.  This  effect  can  be  attributed  to  the  fact  that  the 
sound  velocity  in  the  atmosphere  is  temperature-dependent, 
and  also  attributed  to  the  effects  of  atmospheric  winds. 

Thus,  the  combined  variations  of  temperature  and  wind  veloc¬ 
ity  with  altitude  give  an  effective  corresponding  variation 
of  sound  velocity  with  altitude.  Refraction  of  sound  nd 
shock  rays,  therefore,  occurs.  Rays  initially  parallel  to 
the  earth's  surface  are  refracted  upward  under  conditions 
where  the  velocity  decreases  with  altitude  and  rays  are  re¬ 
fracted  downward  when  velocity  increases  with  altitude.  In 
this  way,  the  atmosphere  can  act  as  a  lens  on  shock  or  sound 
waves  and,  under  appropriate  conditions,  can  focus  or  concen¬ 
trate  the  shock  energy  on  an  area  with  a  resulting  significant 
increase  in  shock  strength. 

While  it  is  difficult  to  imagine  a  practical  appli¬ 
cation  for  atmospheric  acoustic  focusing,  the  concept  of  an 
inhomogeneous  refractive  medium  has  been  used  to  direct 
acoustic  energy  under  water. **  A  spherical  focusing  device 
was  designed  so  that  the  acoustical  refractive  index  varies 
with  the  radial  distance  from  the  center  in  the  proper  manner 
to  collimate  or  focus  acoustic  rays  entering  the  device.  A 
potential  problem  associated  with  this  approach  in  air  is 
that  some  means  to  contain  the  media  with  varying  refractive 
index  must  be  used  that  would  have  a  high  density  relative  to 
air.  The  resulting  acoustic  impedance  mismatch  could  cause 
a  high  reflection  coefficient  for  the  wave  as  it  impinged  on 
the  lens,  with  a  corresponding  transmission  loss  in  passing 
through  the  lens.  The  concept  of  impedance  for  shock  waves 
is  relatively  new  and  apparently  no  investigations  of  the 
feasibility  of  employing  the  Luneberg  lens  for  focusing  shock 
waves  in  air  have  been  reported.5 
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C.  Propagation  of  Planar  Shock  Fronts 

Another  question  that  must  be  considered  relative 
to  the  feasibility  of  directing  shock  waves  is  the  following: 

If  a  planar  shock  of  small  finite  area  could  be  formed,  over 
what  effective  distances  might  it  be  expected  to  propagate 
intact?  Some  degradation  losses  at  the  edges  of  the  shock 
as  it  propagated  would  certainly  be  expected. 

To  gain  some  insight  into  the  extent  of  such  losses 
on  a  propagating  planar  shock  wave,  a  hypothetical  problem 
was  formulated  and  Whitham's  theory  was  used  to  investigate 
the  subsequent  shock  behavior. 

The  two-dimensional  version  of  Whitham's  theory 
associates  successive  positions  of  a  curved  shock  moving 
through  a  medium  by  a  set  of  curves .  The  orthogonal  tra¬ 
jectories  of  these  curves  are  called  rays  and  are  analagous 
to  the  rays  of  geometrical  acoustics.  The  set  of  curves  and 
trajectories  corresponding  to  the  movement  of  a  shock  through 
a  medium  are  used  as  the  basis  for  an' orthogonal  curvilinear 
coordinate  system  so  that  the  shock  positions  and  rays  corres¬ 
pond  to  curves  of  constant  a  and  constant  B,  respectively 
(Figure  3).  The  basis  of  the  therry  is  that  a  shock  is 
always  perpendicular  to  its  direction  of  propagation  as  it 
moves  in  a  uniform  atmosphere  at  rest  and  that  the  localized 
shock  Mach  number  is  related  to  the  area  of  the  shock  front 
"trapped"  between  two  rays  of  the  shock. 

The  coordinate  o  is  conventionally  taken  to  be 
o  ■  at,  where  a  is  the  sonic  velocity  in  the  undisturbed  gas 
ahead°of  the  shock.  The  distance  along  a  ray  between  shocks 
located  at  a  and  da  is  then  M(a,B)da  where  M  is  the  Mach 
number  at  a,B.  Similarly,  the  corresponding  spacing  between 
rays  at  B,  B+dB,  is  A(a,B)dB°  It  can  be  shown  that  geometrical 
considerations  alone  require 


that 


36  m  1  3A 
7F  M  Ta 


and 


30  1  3M 

7a  IT  7F 


(1) 


(2) 
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where  0(o,g)  is  the  angle  made  by  the  ray  at  a,  £5.  The  angle  0 
can  be  eliminated  from  the  above-cited  two  equations  to  give 
a  differential  relation  between  M  and  A.  An  additional  re¬ 
lation  is  obtained  from  consideration  of  the  energy  flux  down 
a  ray  channel.  The  energy  flux  is  proportional  to  the  square 
of  the  pressure  amplitude  multiplied  by  the  channel  area. 
Within  a  given  channel,  the  flux  is  constant  so  that  the 
strength  is  proportional  to  For  weak  shocks,  the  shock 

strength  is  proportional  to  M-l,  so  one  has 

M-l  «  A_J*  (3) 


as  the  second  relation  between  M  and  A  that,  in  principle, 
allows  solution  for  M  as  a  function  of  a  and  0. 

In  adapting  the  problem  to  the  computer,  Equation  (2) 
was  used  in  the  computation  of  each  successive  new  shock 
position,  while  Equation  (1)  was  satisfied  indirectly  by 
recalculation  of  the  ray  direction  normal  to  the  new  shock 
position  at  each  step.  The  relation 

M-l  =  A*  **  (4) 


was  used  to  calculate  the  Mach  number  corresponding  to  the 
area  between  adjacent  rays  at  each  new  shock  position.  A 
typical  FORTRAN  IV  program  used  is  listed  in  the  Appendix. 

The  planar  problem  investigated  involved  a  shock 
wave  propagating  out  from  a  wall  (Figure  4).  The  problem  is 
hypothetical  in  the  sense  that  nothing  is  said  concerning  the 
source  of  energy  driving  the  shock.  The  shock  has  radial 
symmetry  and  consists  of  a  central  planar  section  that  is 
five  feet  in  radius,  surrounded  by  a  curved  section  that  forms 
a  boundary  with  the  wall.  The  curved  section  serves  to  facil¬ 
itate  the  analysis  by  the  avoidance  of  a  discontinuity  at  the 
edge  of  the  shock.  While  it  may  be  argued  that  the  shock  is 
not  truly  planar,  the  degradation  losses  of  the  shock  con¬ 
sidered  should  be  less  severe  than  those  of  a  "true  planar 
shock"  by  virtue  of  the  sharp  pressure  discontinuity  present 
at  the  edge  of  the  latter. 

A  plot  of  shock  Mach  number  at  the  center  line  as  a 
function  of  distance  from  the  initial  shock  position  for 
three  different  initial  shock  sizes  is  illustrated  in  Figure  5. 
As  can  be  seen  from  the  figure,  the  center  line  Mach  number 
remains  constant  at  the  value  set  in  as  an  initial  condition 


9 


10 


FIGURE  4  INITIAL  H 
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until  the  wave  has  propagated  a  distance  approximately  80 
per  cent  of  the  initial  planar  section  diameter.  Beyond  this 
point,  the  center  line  Mach  number  decays  quite  rapidly  with 
a  shape  that,  at  least  for  the  smallest  shock,  strongly  re¬ 
sembles  the  decay  of  a  spherical  blast  wave. 

The  behavior  of  the  shock  profile  as  a  whole  is 
depicted  in  Figure  6.  The  tendency  of  the  "planar"  shock 
toward  a  progressively  more  spherical  shape  is  quite  evident. 

II.  Concave  Shock  Fronts 

A.  Analysis  Based  on  Whithorn' s  Theory 

The  general  type  of  behavior  to  be  expected  of  a 
concavity  in  a  shock  front,  based  on  Whitham's  theory,  is 
indicated  in  Figure  3.  As  this  figure  indicates,  formation 
of  a  caustic  does  not  occur.  Instead,  as  the  front  propa¬ 
gates,  the  rays  begin  to  converge  with  a  resultant  strengthen¬ 
ing  of  the  front  in  the  region  of  the  concavity.  The  strength¬ 
ening  in  this  region  increases  the  velocity  there  and  tends 
to  smooth  out  the  shape  so  that  any  intersection  of  the  rays 
does  not  occur. 

Although  focusing  in  the  usual  sense  of  geometric 
acoustics  does  not  occur,  the  concavity  does  subsequently 
create  a  region  of  increased  "ray  density"  that  implies  a 
region  of  increased  shock  strength.  To  investigate  the  po¬ 
tential  ranges  and  relative  strengths  that  could  be  obtained 
by  virtue  of  this  effect,  several  shock  waves  were  analyzed 
having  various  types  and  degrees  of  concavity.  The  basic 
method  used  was  similar  to  that  used  in  the  analysis  of  planar 
shocks . 


The  effects  of  varying  degrees  of  spherically  shaped 
concavity  on  the  center  line  Mach  number  as  a  function  of 
distance  from  the  initial  shock  position  are  illustrated  in 
Figure  7.  Corresponding  center  line  behavior  for  an  initially 
planar  shock  of  the  same  initial  size  is  included  for  reference 
purposes.  The  effects  of  increased  radius  of  curvature  in  the 
concave  section  are  to  lower  the  peak  center  line  Mach  number 
attained  and  to  increase  the  distance  at  which  the  peak  occurs. 
The  planar  curve  apparently  can  be  thought  of  as  representing 
the  limiting  case  of  infinite  radius  of  curvature. 

The  effects  of  elliptical,  parabolic,  and  hyperbolic 
concave  curvatures,  respectively,  on  the  center  line  Mach 
numbers  as  a  function  of  distance  are  shown  in  Figures  8.  9, 
and  10.  As  in  the  spherical  case,  the  general  effect  of  in¬ 
creased  curvature  was  to  increase  the  degree  to  which  the 
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FIGURE  7  CENTER  LINE  RESULTS  FOR  PROPAGATION 
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-FIGURE  8  CENTER  LINE  RESULTS 

FOR  PROPAGATION  OF  INITIALLY  ELLIPTIC  CONCAVE  SHOCK 
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DISTANCE  TRAVELED  IN  UNITS  OF  SHOCK  DIAMETER 

FIGURE  10  CENTER  LINE  RESULTS  FOR  PROPAGATION 

OF  INITIALLY  HYPERBOLIC  SHOCK 


Mach  number  is  peaked  and  to  shorten  the  distance  from  the 
initial  position  at  which  the  peak  occurs.  The  Mach  peaks 
seemed  most  pronounced  for  the  elliptically  shaped  initial 
configurations.  In  all  cases,  however,  no  significant  gain 
in  Mach  number  over  that  of  the  initially  planar  configur¬ 
ation  was  observed  at  distances  in  excess  of  a  few  diameters 
of  the  initial  shock. 

B.  Experimental  Data  on  Weak  Shocks 

It  was  mentioned  earlier  that  recent  data  reported 
by  Beasley,  et.al.,  relative  to  the  focusing  of  N-waves  has 
indicated  that  very  weak  shocks  can  be  focused  according  to 
the  principles  of  geometrical  acoustics.  The  acoustics  solu¬ 
tion  is  based  on  the  premise  that  rays  normal  to  a  wave  front 
remain  fixed  in  direction.  Consequently,  rays  normal  to  a 
concavity  eventually  form  a  caustic,  or  focus.  Beasley's 
group  investigated  the  focusing  behavior  of  N-waves  having 
amplitudes  between  1.48  and  4.61  pounds  per  square  foot  and 
having  wavelengths  short  in  comparison  to  the  dimensions  of 
the  focusing  mirrors  used.  Their  results  indicated  that, 
when  such  shocks  were  reflected  from  a  focusing  mirror,  they 
passed  through  a  line  or  point  focus  for  the  respective  two 
and  three  dimensional  cases,  with  an  accompanying  phase  re¬ 
versal.  A  sequence  of  Schlieren  photographs  showing  suc¬ 
cessive  stages  in  the  process  was  included. 

The  theory  of  geometrical  acoustics  and  Whitham's 
theory  of  shock  are  basically  quite  similar  in  concept.  How¬ 
ever,  the  effect  of  the  basic  linear  premise  of  geometrical 
acoustics  is  to  represent  all  points  on  a  wave  front  as  having 
the  same  propagation  velocity.  This  leads  to  the  marked  dif¬ 
ference  in  the  focusing  behavior  predicted  by  the  two  theories. 
The  fact  uhat  very  weak  shocks  apparently  conform  to  the 
principles  of  geometric  acoustics  is  quite  significant.  Yet, 
at  higher  amplitudes,  the  increased  deviation  of  the  shock 
velocity  from  the  implicitly  assumed  sonic  velocity  of  geo¬ 
metric  acoustics  must  begin  to  become  apparent  in  focusing. 

The  amplitudes  at  which  these  changes  in  focusing  behavior 
occur  and  the  nature  of  any  transitional  focusing  phenomena 
would  be  a  significant  contribution  to  the  concept  of  shock 
focusing. 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  propagation  of  planar  and  concave  shock  wave'  sec¬ 
tions  was  studied  on  the  basis  of  the  shock  theory  of  G.  B. 
Whitham.  The  results  indicate  that  the  feasibility  of 
directing  shock  waves  of  small  frontal  area  over  significant 
distances  without  spreading  losses  is  very  much  in  doubt. 
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However,  experimental  data  on  the  focusing  of  very  weak  shocks 
indicate  that  they  do  not  conform  to  Whitham's  theory.  Con¬ 
sequently,  it  is  recommended  that  the  apparent  discrepancy  in 
the  theory  of  weak  shock  focusing  behavior  be  resolved  ex¬ 
perimentally  before  a  final  conclusion  is  reached  on  the 
feasibility  of  focusing  this  type  of  shock  wave. 
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APPENDIX 


SHOCK  WAVE  PROPAGATION  PROGRAM 
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JOB  I  »,» 


//-oiMf.H  JOB  (01740  39,  3,3000  ), 'AW COR -REMOTE- 02-SAVE' 

(. 

0  NOMENCLATURE 

r. 

(.  K  A  =  NUMRFR  OF  FORWARD  STEPS  FOR  SHOCK 

(  X (  1  )  =  AX  1 AL  DISTANCE 

f.  R(  I  )  RAOIAL  DISTANCE 

C  M (  I  )  =  MACH  NUMBER  ON  CHORD  L  INF  CONNECTING  POINTS  I  f.  1  +  1 

(.  A  (  I  )  AREA  OF  CHORD  FROSTRDM  /PI 

r  DB( I )  =  LENGTH  OF  CHORD  LINE 

(  DP ( I )  =  PRESSURE  RISE  ASSOCIATED  WITH  M ( I ) 

r.  HS  =  SOM  OF  DB(  I  )  S 

(.  AVR  =  RS/30,  THAT  IS,  INCREMENT  LENGTH  TO  CALCULATE  NEXT  FORWARD  STEP 

C  TH  =  ANGLE  OF  SHOCK  RAY  TO  +X  DIRECTION 

(.  OMOR  =  RATE  OF  CHANGE  OF  MACH  NUMBER  WITH  DISTANCE  ALONG  CHORD 

C  MAV  =  AVERAGE  MACH  NUMBER  AT  LATTICE  POINT 

C  THT  =  ACTUAL  DIRECTION  USED  TO  CALCULATE  LOCATION  OF  NEW  POINT, 

(.  =  THFTA  -t-DELTA  THETA/? 

C  MA  =  2.0  DIVIDED  BY  MACH  (1)  TO  NORMALIZE  AVR 

(. 

t. 

IMPLICIT  RFAL*R( A-h.m-7 ) 

DOUBLE  PRECIS'ON  UATAN,DSORT,DCUS, OSIN, DABS 

DIMENSION  X<  3 .  ) ,R( 31  )  ,M( 30) ,A( 30 ) ,DR( 30) ,DP( 30) ,Y<  31) ,S( 31 ) 

KA  =  400 

HO  ?H0  K  =  l  ,  K  A 

C  IF  K  =  1  SET  INITIAL  VALUES  OF  X , R , AND  M 
I  F  <  K  —  l  >10,10,140 

r.  SET  VAI.OFS  of  X,  R  AND  M  FOR  I  =  27  THRU  30 

10  DO  pn  I  =  1,4 

Al=?«I— 1 

R  (  1  +2  4  )=S.DO+O.BOO*DS IU( 3.14159U0*AI/20.D0  ) 

X ( 2  7  )  =  4.9«3M44?D0 

X ( 2  H  )  =  4.943303400 

X(?9)  =  4. R 3 33 3 3300 

X ( 30 )  =  4. 72(3993700 

20  i-(I+?M  =  2,00  -  0,07777DD*AI 
f.  SET  VALUES  OF  X  ANO  M  FOR  I  =  1  THRU  23 

OO  30  I  =  1,2ft 
X  (  I  )  =  3.000 

30  3(|)  =  2.00 

S  A R C  =  3 . 1 M  3900/33. DO 

(  SET  VALUES  OF  R  FOR  I  =  1  THRU  2(3 

DO  /VO  I  =  1,4 

A  I  =  ?« I 

40  R  (  T  +  2  2  )  =  3.00  +  SAVC  »  ( A  I  -  9. DO) 

'  SMOOTH  OUT  R  FDR  I  =  2  THRU  ?? 

I  S=o 

ID  30  I  =  1,22 

3l '  I  >1  =  (9  +  1 

'13=19 

I 
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AI)FLT*(R(23)-44.I)0*SARC>/AIS 

r<  i )  =  o.no 
on  (so  i  *  1,21 

A  1  =  I 

40  R  ( 2  '4  -  I)  =  R<24  -  I)  -  2.00  *  SAHC  -  AllELT*AI 
MP2S=( 3.00*M(24)+M<2H  )  )/4.00 
MP27=<  3.D0*M(  2R)*M«  24)  )/4.l)f) 

M(2S)=(M(2S)  +MP2S  )  /  2  •  DO 
M(  2  7  )  =  <M(  2  7  )+MP27  )/2  .1)0 
MP2A=(M(24)+rt<2«)  )/2.l)0 
M<24)s(M(?4)+MP24)/2.l>0 
r.  SFT  FINAL  VAI.OF S  OF  «  AMO  X  FOR  I  =  31 

X  (  31  )  =  4.  SI  10 
R  (  31  )  *S  »  SOO 

0  HFSFT  X  VAI.IO-S  FOR  SPECIAL  SHARKS 

70  1)0  HO  1  =  1  ,2(S 

X  I  1  )  sOAM  S  I  I  I0.7SI)OMOSOHT<  2S.D0-R  (  1  )#*  2  )  l/S.OO)  -S.OO) 
ho  r, timt  i  nof 

r,si  .400 
WRITF  14,40) 
or.  FORMAT (MHOX  a  1) 

•-RITF  14,100) 

loo  FORMAT ( •  X  R  MACH. MO.  OFLTA  PR 

1FSSORF  OH*//) 

r.  CALCOLATF  All),  OKI  I),  AMO  OR  (  I  ) 

on  no  r  =  1,30 

A(l)  =  (R(I  +  l)-fR(i)  )«OSORT(  (  X  (  I  )-X  (  1  *-1  )  )*«2+(K(  I  )  -R  (  ]  +  l  ))*#?) 

OH ( 1  )=OSORT ( <X< I )-X< 1  +  1 )  )**?+<R< I  )-K ( 1+)  )  )**2  ) 

or(  I  )=2.oo»r;#  (  m<i  )m#2-1  .no)  /  (t;+i  .oo  ) 

110  WRITF  (4, 120)1, X(  I  ),K(  1  ),.>((  I)  ,  OR  (  I  )  ,  OH  (  I  ) 

120  FORMAT ( I  4, SOI S.H ) 

"Kl  TF  <  4 , ) 30 ) X ( 3 1  ) , R  (  3  1  ) 

130  FORM"  T  ( 4H  3 1 , 2 0 1 S . H ) 

’-'A  =  (  7.00/  (  S.OOWMI  1  )  )  ) 

(01  TO  2H0 

0  r.AI  r.Ol.ATF  FORWARD  I  NCKFMPmI  fM  ALPHA 
140  US  =  0.00 

Oil  ISO  1  =  1,30 
ISO  MS  =  l<S  +  OH  (  I  ) 

I F (K. LF.200  )  AVH  =  (HS/3O0.D0)  »  MA 

I  F  (  K  .0T.200.AVD.K.LF  .SO()  )  AVH  =  (HS/ISO.OO)  a  MA 
I F (K.OI .SOO. ANO.K.LF.l 000 ) AVH  =  (HS/IOO.DO)  *  MA 
I F ( K .fiT . 1000. ANO.K.LF .1S00  )  AVH  =  (HS/SO.DO)  *  MA 
I  F  (K.liT.  1  SOO  )  AVH  =  (HS/30.00)  «  mA 
C  CALCOLATF  LOCATIOM  OF  MFW  PoJMTS 

f.  RFRI.ACF  X  AMO  R  VAI.OFS  WITH  Y  AMI)  S  OF  S  I  ON  A  T  IONS 

140  I  HI  170  I  =  1,31 
Y (  1  )  =X (  I  ) 

sm=«  m 

170  COMTIMIIF 

C  CALCOLATF  I H ,  OMOH,  MAV,  1HI,  AMO  XII)  AMO  R(I) 

1  HO  1)0  200  I  =  2,30 

TH=II)ATAM(  ( Y(  !-l  )-YI  1  )  )/(S( I  )-SI  1-1  )  )  )+DATAN<  ( Y (  I ) -Y (  1  +  1)1 
l/ISI 1+1 )-S ( 1)1) )/2.DO 
0M0H  =  2 .00#  (  M(  I  ) -h (  I  -1  )  )  /  (  |)H  (  I  )  +0K  (  !  - 1  )  ) 

MAV=(M(  |  )*0H(  I-H+lXII-1  )*I)H(  I  )  )  /  (  OH  |  I  )  +0H  (  I  -  l  )  ) 

THT  =  TO-AVH»MAV#OMOH/2.I)() 

1  VO  X  <  I  )  =  XII)  +  AVH  #  MA V  *  DCOSITHT) 

R  (  I  )  =  R  (  I  )  +  AVH  *  mAV  #  OSIM1HT) 

200  COMTIMIIF 
X  (  1 ) *X ( 2  1 
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K  (  HI  )  =R  (  30  )  +0.0011)0 

(.  R  (  31  )  SFT  LARGFR  THAN  R  (  30  )  Til  STOP  A7AN  ARGUMENT  FROM  =  INFINITY 

(.  CALCULATE  AT,  M(I),  AI,  l)H  (  I  )  ,  AND  OP  (  I  ) 

2  l  o  on  ?  ho  i  =  1,30 

A  I  =  (R  (  1  +1  )+R(  I  )  )*l)SORT(  (  XII  l-XII  +  l  )  )#*2+<R  (  I  )-R  (  1  +  1  )  >#*2  ) 

??{)  M (  I  )  =  1  .r>0+(  M(  I  )-l  .1)0  )  *  OSCJRT  (  A  (  I  )  / AT  ) 

A ( I ) = AT 

I  )R|  I  ) sOSURT (  ( X(l)-X(  1  +  1 ))**2  +  <K  (  I  )-K ( 1+1  )  )**2  ) 

2  30  i)P(t)  =  2  .  I)0*G*  (Mil)  +#2  -  1 • DO  )  /  (  G  +  1.1)0) 

C.  WR1TF  (HIT  THF  WAVE  FRONT  POSITION  AND  PROPERTIES 
IF  IK. IF. 299)  Gil  TO  2H0 
232  I F( IK/ 10 )*10-K) 280,240,280 
2M1  WHITT  I  ft , 2  50  )  K 
250  FORMAT  I  7 HO  K  =  ,IH) 

WRITF  (ft, 100) 

DO  2ftO  I  =  1,30 

?ftO  WRITF  (ft,  120)  I,  XI  I  )  ,  R  I  I  )  ,  M  (  I  ),DP(  1  )  ,  IJB  I  I  ) 

WRITF  (ft, 1301X131 ) ,R I  3 1 ) 

GO  TO  2 HO 
2R0  CONTI  WMF 
310  STOP 
FMO 

5  F  M  T  R  Y 
5  STOP 
/* 

II 
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